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^vq We are doubtlessly familiar with some edition of Jackson's tome on electrodynamics, and 

Schwinger's calculation of the anomalous magnetic moment of the electron in QED. From the 
perspective of strong interactions, however, electromagnetic effects usually amount to negligible 



contributions. Despite this fact, electromagnetic probes have always been a fundamental source 
for our knowledge of QCD experimentally. Elastic scattering of electrons off nucleons provides 
us a window to their distributions of charge and magnetism. To account for the spectrum of QCD 
t-H at the percent level, moreover, we need isospin breaking introduced from both quark masses and 

electric charges. This overview concerns some of the prospects and progress of studying electro- 
magnetic effects in QCD. Our focus is divided between classical and quantum effects. In classical 
electromagnetic fields, the dynamical response of QCD to external conditions can be investigated. 
The vacuum and hadrons alike should be viewed as media which respond to external fields: both 
magnetize and polarize in magnetic fields, for example. At the quantum level, electromagnetism 
and QCD renormalize each other. In the era of high precision lattice computations, both strong 
and electromagnetic contributions must be accounted for to make predictions at the percent level. 
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1. Scope 

Electrodynamics in both its classical and quantum form is a staple in the study of physics. 
While the electromagnetic interactions of quarks are dwarfed by their strong interactions, electro- 
magnetic observables and small electrodynamic effects can ultimately give one an intuitive picture 
of the underlying QCD environment present in the vacuum and in hadrons. In the era of high pre- 
cision lattice QCD, moreover, precise quantitative understanding of the strong interaction requires 
accounting for electromagnetic effects. The purpose of this overview is to touch briefly upon a 
few topics involving QCD in the presence of electromagnetism. In the past few years, the level of 
activity in this field has grown considerably. Various different electromagnetic effects are currently 
being addressed using lattice QCD techniques; and, we expect many of these calculations to be 
refined in the foreseeable future. We do not attempt a complete review of this subject. Our goal 
is to give one a general picture of the prospects and progress in the field by highlighting just a 
few areas. In many of these areas, a detailed technical review of the subject is not yet warranted. 
We look forward to the many future developments that will likely necessitate a thorough review of 
this field. Our overview is divided into two parts: one dealing with aspects of QCD in classical 
electromagnetic fields, and a second dealing with aspects of QCD in the presence of QED. 

2. Lattice QCD with Classical Electrodynamics 

Studying QCD in the presence of classical electromagnetic fields allows one to address how 
QCD responds to external conditions. We will focus on just two aspects of this response. First, we 
look at how hadrons respond to weak electric fields. Second, we consider how the vacuum of QCD 
responds to strong magnetic fields. 

2.1 Weak Electric Fields 

Materials are characterized by two broad types of properties. Materials have intrinsic prop- 
erties that are determined by their underlying structure. On the other hand, materials also have 
extrinsic properties that are determined by their response to external conditions. For hadrons, very 
little is known about their extrinsic properties. The only extrinsic properties reported in the PDG 
concern the response of the nucleon to external electromagnetic fields. This response is captured 
by the electromagnetic polarizabilities. 

Electric polarizabilities are familiar from classical electrodynamics. A material without an 
electric dipole moment will nevertheless polarize in an applied electric field. The induced dipole 
moment, d-E, will be proportional to the applied field and oppositely oriented, d£ = —(XeE, with a 
coefficient, Cfe, representing the electric polarizability of the material. Consequently the effective 
Hamiltonian for the material in an external field acquires a term at second order in the electric field, 
AH = —j(XeE 2 , with higher-order terms assumed to be negligibly small. Electric polarizabilities 
of atoms and molecules play a crucial role in chemistry and biology. Quantum mechanically one 
can compute the polarizabilitiy of the ground-state hydrogen atom exactly, a E I ~ atom = |^ (^7la 3 B ) °c 
e~ 6 m~ 3 . Not surprisingly the polarizability scales with the volume of the hydrogen atom ground 
state with a coefficient that is order unity. In the limit of vanishingly small electric charge, the 
hydrogen atom becomes extremely loosely bound, and the electron cloud is readily polarized by 
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the applied field. Furthermore, in the limit of vanishingly small electron mass, the hydrogen atom 
becomes completely diffuse. 

For hadrons, we thus have considerable intuition about their polarizabilities, the form of which 
must scale with the typical hadron's volume: a^ adro " = jV(Xf. s . (^7T [fm 3 ]), with ,jV a hadron- 
dependent pure number. In the limit of vanishingly small electric charge, the electric field decouples 
from the quarks. The scale of the hadron's volume, however, is a dynamically generated quantity 
depending on non-perturbative QCD. It is consequently less obvious how the electric polarizability 
should behave near the chiral limit, m q — > 0. In the MIT bag model [1], the nucleon polarizability is 
insensitive to the quark mass, Cfe ~ e2 ^iag m tne cri iral limit; whereas, in chiral perturbation theory, 

i —1/2 —1/2 _? 

the polarizability has a striking non-analytic behavior, Cfe ~ e m q ^q CD A x , associated with the 
charged pion cloud [2]. In the chiral limit, hadrons should just barely become diffuse; but, whether 
the physical values of the light quark masses put us in the regime where this term dominates is 
unclear. For the charged pion, predictions from chiral perturbation theory for polarizabilities are 
almost a factor of two smaller than the most recent experimental determination [3]. 

Electromagnetic polarizabilities are accessible in Compton scattering experiments, and there 
are on- going experimental efforts to determine these quantities precisely. The COMPASS experi- 
ment at CERN hopes to provide the final word on charged pion polarizabilities, and the first word on 
charged kaon polarizabilities. Low-energy Compton scattering on deuterium has been performed 
recently at MAX-Lab in Lund to allow for the extraction of neutron polarizabilities. A significant 
effort is underway using HlyS at TUNL to measure Compton scattering off hydrogen and deu- 
terium targets. It is thus timely for the lattice to compute electromagnetic polarizabilities to help 
resolve long-standing discrepancies. The ability to vary the quark masses is a virtue, moreover, as 
it allows for a crucial test of low-energy QCD from first principles. 

The Compton scattering tensor between hadrons is a four-point function, namely T^ v (k,k') = 
J x e' k ' x ~ y (H\^ {j^(x)J v (y)} \H), which lattice techniques have not yet tackled. Even if the 
full four-point function could be evaluated on the lattice, access to polarizabilities would only be 
through long extrapolations to vanishing photon momentum. Fortunately there exists an alternative 
that dates to the early period of quenched lattice QCD: the background field method. Background 
magnetic fields were used in the first calculations of nucleon magnetic moments [4], and shortly 
after it was realized that effects at second order in background electric fields could be used to 
extract electric polarizabilities of neutral hadrons [5]. 

The background field method was revisited a few years ago to determine magnetic moments 
and magnetic polarizabilities of hadrons, and electric polarizabilities of neutral hadrons [6]; these 
studies were completely quenched. Non-uniform external fields were employed by these studies 
with effects from field gradients mitigated by imposing Dirichlet boundary conditions in space 
(space and time) for magnetic (electric) fields. The corresponding eigenstates are standing waves, 
which were not addressed in these works. The boundary conditions, moreover, introduce non- 
perturbative effects that are difficult to quantify. 

Further computations have utilized external electric fields to deduce the electric dipole moment 
of the nucleon in the presence of the (9-term [7]. These computations are partially quenched (in 
this case, they are quenched only by shutting off the electric charges of the sea quarks), but still 
suffer a spike in the electric field at the temporal boundary. Here standard antiperiodic boundary 
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Figure 1: Effective mass plots for (connected) neutral pion correlation functions in external fields. On 
the left, we compare zero field and the naive implementation of an electric field. On the right, we compare 
quantized, n = 3, and non-quantized, n = e = 2.71828 . . ., values of the electric field with a boundary twist 
included. In both cases, the correlators have been translated in time so that t — 12 corresponds to the edge 
of the lattice. Plots are taken from [13]. 



conditions are imposed on the quarks in time, with sources and sinks located as far as possible 
from the spike. Additional studies have given preliminary reports on the pion mass dependence 
and volume dependence of neutral hadron electric polarizabilities, and a comparison of boundary 
conditions [8]. 1 We summarize three recent developments in external electric field computations. 

1). To apply a uniform electric or magnetic field on a torus, quantization conditions must be 
met [10]. The generalization for a discrete torus, i.e. a lattice, has been detailed in [11]. For 
example, the allowed strengths of uniform (Euclidean) electric fields $ satisfy 




(2.1) 



with q as the quark electric charge, and n an integer. 2 Lattices are currently large enough to support 
perturbatively small electric fields. Effects of quantized versus non-quantized fields were studied 
in [13]. For a neutral pion in a uniform electric field, the energy has the form, E n o = m n o + 
iOgV 2 , up to higher-order corrections in the field strength. With a non-uniform field, there are 
a tower of additional terms involving gradients. Each term in the tower is accompanied by an 
unknown coefficient. The electric field $z, can be specified by a time-dependent vector potential 
= (0,0, — <fx4,0). To include this field on the lattice, one uses Abelian gauge links, Up = 
exp(i^A jU ). This naive implementation leads to a large field gradient between X4 = J3 — 1 and*4 = 0. 

'it is of considerable interest to note an alternate approach that is complimentary to the background field method. 
Instead of including background fields in the action and later hunting for effects in two-point functions that are a desired 
order in the strength of the external field, one could expand correlation functions from the outset to isolate physics at the 
desired order. This alternate approach has been developed and demonstrated by measuring the electric polarizability of 
the neutron including all disconnected terms for m K = 760 MeV [9]. A salient feature of this approach is the ability to 
treat explicitly volume effects from the holonomy of the external field. A downside is that non-uniform fields are used 
to perturbatively expand correlators. 

2 While quantized field strengths eliminate the gradient systematic, there are still non-trivial finite size effects. These 
arise due to non-vanishing Polyakov loops, e.g. the naive electric field has the holonomy <I>3 = Jq dxiAs(x) = —SLx^. 
At the level of hadrons, charged pions propagating around the world /7-times must be dressed with a (spatial) Polyakov 
loop, e'"* 3 . Finite volume effects thus introduce a time-dependent potential, however, coefficients in this potential are 
exponentially suppressed, « e -nm K L f or aS y m pt tically large volumes, see [12]. 
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Figure 2: Effective mass plot for the charged pion correlator calculated in an external electric field, with 
strength n = 3. A fit to the correlator is shown allowing for excited-state contamination, and the extracted 
"rest energy" is plotted as the solid red line. Plot taken from [16]. 

A very dramatic effect can be seen on the (connected) neutral pion correlation function, as shown in 
Figure 1. The field gradient can be mitigated by including a boundary twist, [7' WMt = exp^qA'^'), 
with A , ^' ,st = (0,0,0,j3(fx35 :r4 j3_i). The resulting field is uniform except for a huge gradient at the 
corner of the lattice, {xt,,x 4 ) = (L — I,j8 — 1). This gradient can be eliminated when the field is 
quantized as in Eq. (2.1), see Figure 1. 

2). Charged particle properties can be deduced in external electric fields. This was first sug- 
gested in [14], where a non-relativistic method was developed. In light of the size of the external 
field strengths, a relativistic generalization is needed in practice, and can be found in [15]. The 
basic idea behind the treatment of charged particle correlators, is to write down an effective action 
for the hadron in the external field. 3 For example, a composite scalar of charge Q in an external 
electric field, = (0,0, — £x 4 ,Qi), has an effective action, 



S£ e ff (p = 0,x 4 ) = <t>Hx 4 ) 



(j>{x 4 ), (2.2) 



where the "rest energy," E($) = m§ + ^a^S 12 -\ , includes non-minimal couplings such as the 

polarizability. The propagator for Eq. (2.2) can be expressed easily using Schwinger's proper-time 
trick, 

1 1 p 00 

Gto = ( 1^, eJ O) = o / dse-**(<u\e-"\0). (2.3) 



2Jf + £ 2 ' ' 2 Jo 

Charged pions and kaons were investigated in external electric fields [16], and their two-point 
correlators were fit using Eq. (2.3). Extracted values for E($ ) at various field strengths give one the 
ability to determine charged particle electric polarizabilities. An effective mass plot for a charged 
pion correlation function is shown in Figure 2. 

3). Spin-i particles have also been treated in electric fields. There is a subtlety even for the 
neutron because of the spin interaction with the external field. Using the neutron effective action, 
££ e ff = [y^d^ +E{S) — ^cfyvF^v] N, one finds that unpolarized neutron two-point functions 



3 There is also implicitly the assumption that the spectrum of lowest-lying states remains ordered as in zero field. 
This assumption can break down in considerably strong fields, e.g. S ~ Mjr, where M* is an excited state mass. 
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are given by the standard exponential, Zexp(— E e f/z), with a falloff determined by E e ff = M + 
\S 2 [a,E — 1 ■ The additional term is a ~ 10% effect at the physical pion mass, but we should 
make no assumption of the effect at larger-than-physical values of the pion mass. If we can find 
a way to remove this term, moreover, we will have succeeded in determining both the electric 
polarizability and magnetic moment of the neutron. This problem was discovered in [17], where 
a solution was posed and verified using lattice QCD. By analogy with a magnetic field, where 
5C7m V Fm V = S B and one forms spin-projected correlation functions to access magnetic moments, 
in an electric field, we have h<7^, v Fu V = K E, where the matrices K are the spin-^ generators of 
boosts. In forming boost-projected correlators, i.e. by tracing with ^(1 ±£3), one finds G±(t) = 
Z M ± ^j^J exp(— .Eg//T). Thus a simultaneous fit to these correlation functions will lead to a 
determination of jU and o^. The additional magnetic moment squared term leads to a 50% effect at 
m n = 390 Me V. A generalization for charged spin-| particles was also demonstrated in [17]. 

At this conference: There were a number of developments related to the study of hadrons 
in weak external fields reported. A preliminary look at quenched neutron correlation functions in 
uniform magnetic fields was given [18]. Because the charge and anomalous magnetic moment in- 
teractions do not commute for relativistic spin-j particles, Landau levels must be treated to extract 
the magnetic moments and polarizabilities. As chiral symmetry properties of lattice fermions may 
play an important role in polarizability calculations, issues related to using overlap fermions in 
external fields were discussed [19]. When one considers external fields with slowly varying space- 
time dependence, additional quantities, known as spin-polarizabilities, characterize the hadron's 
response. Such quantities are even more tightly constrained in chiral perturbation theory. Feasi- 
bility of extracting spin polarizabilities using background fields was discussed [20]. A promising 
alternate method to compute the electric spin polarizability of the neutron was unveiled [21]. 

2.2 Strong Magnetic Fields 

Lattices are not yet large enough to support several values of perturbatively small, uniform 
magnetic fields. This is because the 't Hooft quantization condition for magnetic fields gives val- 
ues qB = 2nn/L 2 , which are larger than the allowed electric fields in Eq. (2.1), because generally 
one has j3 fa 2L. The dimensionful parameter determining whether a field strength is perturbative is 
the mass of the lightest charged state that couples to the hadron in question. This state must be the 
pion due to QCD inequalities. Whether the pion is the dominant contribution, as in chiral perturba- 
tion theory, is irrelevant. Thus the expansion is governed by e 2 B 2 jm\ ~ n 2 jl ', for m K L 4, and one 
is certainly restricted to n < 3. Notice that e 2 S 2 jm\ is smaller by a factor of j8 2 /L 2 th 4. Using uni- 
form magnetic fields, there has been initial work on the magnetic moments of spin-| baryons [22], 
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Figure 3: Estimates of strong magnetic fields in the universe given in units of Tesla, and m\. 
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and an investigation of the magnetic-field dependence of the proton's wavefunction [23]. Complete 
studies in weak, uniform magnetic fields await larger lattices. 

Rather interesting effects can arise in strong magnetic fields; and, these fields are not only 
accessible to lattices today, they are in a range relevant for heavy ion collisions. The currents 
generated by two beams of large ions (be they gold or lead) give rise to tremendous magnetic fields 
in the region between the two colliding beams. Estimated physical extremes for magnetic fields are 
summarized in Figure 3. Some elements of the conjectured phase diagram for QCD in an external 
magnetic field are depicted in Figure 4. We focus on three regions in the B-T plane. 

1). For \eB\jm\ < 1 and \eB\/A z % <C 1, one can use chiral perturbation theory with strong- 
field power counting to determine how QCD responds. For simplicity, we consider cold QCD, 
i.e. temperatures such that T <C m n . It is known that a magnetic field catalyzes chiral symme- 
try breaking: a magnetic field preserves configurations with particle-antiparticle pairs of oppo- 
site momenta. In the chiral limit, one has the non-analytic enhancement factor, ( YW) b / (WW) = 
1 + log2|<?B|/A^ + fpjjzls calculable in chiral perturbation theory [24]. Low-energy theorems, 
such as the Gell-Mann-Oakes-Renner relation can be generalized in external fields. A quantity 
of particular interest is the magnetization, %, of the vacuum, which we define by, % = lim^o %b, 
where {y^^ivW)b = ZbF^iv- This parameter is of interest for QCD sum-rule calculations, and can 
be accessed experimentally by measuring chiral-odd parton distributions. Unfortunately it is an 
unknown parameter in chiral perturbation theory. 4 Nonetheless, the lattice can fully explore the 
QCD response to strong magnetic fields in this regime and beyond. A pioneering lattice study has 
measured the response of the chiral condensate in quenched SU(2) gauge theory [26], and later 
in quenched QCD [27]. A further study gave the first treatment in two-flavor QCD [28], where 
fully dynamical lattice were generated with staggered quarks coupled to a uniform magnetic field. 
Unquenching the magnetic field was found to have a substantial effect on the behavior of the chi- 
ral condensate. Finally the magnetization of the vacuum was computed in quenched SU (2) gauge 
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Figure 4: A few elements of the conjectured phase diagram of QCD in the B-T plane. We focus on 
questions in regions marked by question marks. Plot adapted from [25], where further details can be found. 



One can use chiral perturbation theory to demonstrate that Xb / X = (w) B / (w) ■ The constant of proportionality, 
however, is the parameter of interest. 
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theory [29], and quenched QCD [27] using a generalization of the Banks-Casher relation. 

2) . A second region of interest in the phase diagram encompasses the transition to the quark- 
gluon plasma. The chiral condensate melts as the temperature is raised, with chiral symmetry 
eventually restored above some temperature. There is another finite temperature transition from 
a confined hadronic phase to a deconfined plasma phase. In zero magnetic field, these transitions 
occur at basically the same temperature. Turning on a magnetic field, however, it is conceivable 
that these transitions might substantially split [30]. Naively chiral symmetry breaking is enhanced 
in an external magnetic field; consequently we might expect chiral symmetry restoration to be 
postponed. If this were the case, there could be a new phase of deconfined, chirally asymmetric 
matter. This possibility was first explored using lattice QCD in [31]. To find the deconfinement 
transition, typically the Polyakov loop is studied. Thus to assess the effect of a magnetic field 
on deconfinement, it is a necessity to have fully dynamical configurations. As mentioned above, 
these configurations were generated with a two-flavor staggered quark action. This study found a 
slight increase in T c , but no noticeable split in the transitions. The transitions appeared to sharpen 
potentially becoming first order. Only one lattice spacing was employed, however, which was 
rather coarse, a « 0.3 f m. Preliminary results from a more complete study were reported at this 
conference, see below. 

3) . Our final look at strong magnetic fields concerns the possibility for novel QCD phase 
transitions. Quenched lattice simulations show evidence for an insulator-(anisotropic) conductor 
transition in strong magnetic fields at non-zero temperatures [32]. As the temperature is lowered, 
there exists the possibility to form an anisotropic superconductor [33]. This is the QCD analogue 
of the electroweak vacuum instability in strong magnetic fields that leads to condensation of weak- 
bosons [34]. In ordinary superconductors, strong magnetic fields will cause the Meissner effect to 
break down. For the present case, however, past a critical value of the magnetic field, B c , charged 
vectors should condense leading to a superconducting state, see Figure 4. 

Unlike the point-like vector bosons of the weak interaction, vectors in QCD are composite 
particles with non-minimal couplings to the external field. The effective action for a charged rho- 
meson in a uniform external field is 

S£ efS = pi [(-D 2 +M 2 p )S a p +D p D a - in P F ap \ p p . (2.4) 

In the action, we have included an anomalous magnetic moment coupling, fi p , but ignored the 
field-strength dependence of the mass, M p , and omitted the width. A charged pion must be one 
of the decay products of a charged rho, but the decay becomes suppressed in large magnetic fields 
because the lowest Landau level for the charged pion has an energy that increases with increasing 
field strength, E n = \/m^ + \eB\. The lowest Landau level for the charged rho, by contrast, has the 
form, E p = \Jm^ — fl p B, when the rho is polarized opposite the direction of the magnetic field. 

In this simplified picture, the vacuum will become unstable past B c = M p j\i p . In this regime, 
a polarized charged rho condensate will form: (0|p||0) = &(x x B), where & is a particular 
periodic function of the coordinates transverse to the magnetic field [35]. Using quenched SU(2) 
gauge theory, quark condensates with charged-rho quantum numbers were found past a critical 
field \eB c \ ~ GeV 2 [36], which suggests that this picture might be correct for QCD. 

At this conference: Results related to the phase structure of QCD in strong magnetic back- 
grounds were discussed. A brief summary was included in the review of QCD at finite temperature 
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and density [37]. The first study of the chiral condensate in QCD with a magnetic field (described 
above) was reviewed [38]. A preliminary report of an extensive study of QCD in external magnetic 
fields was given [39]. In this study, staggered lattices coupled to magnetic fields were generated 
with 2 + 1 flavors at the physical quark masses. Ensembles with different lattice spacings were gen- 
erated, moreover, to enable taking the continuum limit. There was no splitting of deconfinement 
and chiral transitions observed (up to \eB\ ~ GeV 2 ). Unlike the first study, however, the critical 
temperature was found to decrease, with the transitions remaining crossovers. Data from this study 
show that magnetic catalysis of chiral symmetry breaking is not universal: there is an interplay 
between magnetic field, temperature, and quark mass dependence of the chiral condensate. 

Relating to the chiral magnetic effect [40], there were two talks. In [41], the charge separation 
of low modes was investigated in a magnetic field; while, in [42], a so-called chiral chemical 
potential was employed. Previous lattice results on the effect include [43, 27]. We regret that time 
and now space do not permit us a closer look. 

3. Lattice QCD with Quantum Electrodynamics 

In nature, strong and electromagnetic interactions coexist. Consequently they renormalize 
each other. High precision work in lattice QCD requires treatment of electromagnetic effects. This 
has lead to the study of QED corrections to hadron properties. 5 These studies are reviewed after 
we first describe a few QCD applications that require QED. 

1). Accounting for isospin splittings in the spectrum of QCD is one important problem that 
will be solved by treating QED effects on the lattice. In QCD + QED, SU (2)y is broken by the 
light quark mass difference, and quark electric charges. In the nucleon mass, the leading isospin 
breaking is proportional to — m u , and q 2 u — q 2 d . The well-measured nucleon splitting is only a 
fraction of the nucleon mass, (M n — M p )/Mn ~ +0.1%. The electromagnetic mass of the proton 
can be estimated crudely 6 from the classical energy needed to build up a sphere of charge \e\ of 
radius 1 f m, which is oc a,f, s J\ [f m] ~ +1 MeV. Ignoring magnetic effects, the electromagnetic split- 
ting is on the order of (M„ —M p )qed ~ — lMeV, from which we infer the strong isospin splitting, 
(M„ —M p )qcd ~ +2 MeV. To determine the nucleon electromagnetic mass rigorously, one needs to 
account for the distribution of both charge and magnetism within the nucleon. This information is 
provided by the elastic electromagnetic form factors. There are, however, additional inelastic con- 
tributions, and their theoretical status has been discussed recently, see [45] and references therein. 
The electromagnetic mass splitting can be computed using lattice QCD with QED. 

Our discussion of isospin splitting in the nucleon brings up an important point. Currently our 
knowledge of the strong isospin splitting is limited by the removal of the electromagnetic contri- 
bution. This is a general feature. High precision determination of strong interaction parameters, 

5 There are also cases for which QCD corrections are needed for QED observables. The most noteworthy example 
being the high-precision experimental determination of the anomalous magnetic moment of the muon. Confrontation of 
this result with theoretical calculations has been limited by our knowledge of the QCD corrections. The status of these 
non-perturbative corrections to QED computed in lattice QCD was reviewed at this conference [44]. 

6 While this estimate is crude, it is consistent with chiral symmetry constraints. Treating the quark charges as 
arbitrary, the classical sphere estimate has the dependence, (M„ — M p )q£d oc Q% — 0^, = —3(qf, — q^). 
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most importantly: the masses of up and down quarks, is limited by our understanding of QED 
corrections to the inputs. The mass of the up quark is particularly susceptible to uncertainty arising 
from QED corrections to observables because it is the lightest quark. 7 One way to determine strong 
interaction parameters without additional assumptions is to simulate 1 + 1 + 1-fiavor QCD + QED. 

2) . Within the pseudoscalar meson sector, electromagnetic corrections to meson masses are 
universal in the chiral limit. This is known as Dashen's theorem, and it is thus interesting to 
measure its violation. Violation of the theorem can be parameterized by a quantity ^M^ ashen , 
where AMp ashen = (m|+ - m 2 K0 ) QE D - (»4+ - ^ ) Q £ D - By construction, &(a f . s ) terms cancel 
out of AMp ashen , and the leading non- vanishing term has the behavior &(m s GCf^) away from the 
St/ (3) chiral limit. A comparison of lattice and phenomenological results for this quantity was 
shown in [47], with the agreement decidedly poor. 

3) . Beyond masses and mass splittings, isospin breaking effects can complicate the compari- 
son of other strong interaction observables computed on the lattice to experiment. High precision 
determination of the kaon decay constant, fg , is a good example of such a quantity. Experimentally 
the combination ^ x |K. s |/|Kd| is determined to an accuracy of 0.2%. Lattice results on fa/ fit are 
known to 0.4%, using the world average, and are used to extract the CKM matrix elements. At this 
level of precision, isospin breaking effects need to be considered. The experiment measures the 
charged kaon decay width, thus the relevant quantity is f K ± . There are electromagnetic corrections 
to the decay constants of &{otf. s ) which have been accounted for in the analysis used to extract 
CKM matrix elements. The up-quark mass dependence was an additional issue considered in [48]. 
Isospin degenerate simulations extract where m = \(m u +rrid) and the functional depen- 
dence indicates only the lightest valence quark mass, whereas one needs fx (m u ) for the charged 
kaon. Using chiral perturbation theory to account for the valence quark mass dependence leads to 
a ~ 0.2% difference in f K ±. Given the precision of the world average, it is of crucial importance to 
have a high precision determination of m u to account for such effects. One also needs simulations 
at the physical value of m u to end the reliance on chiral perturbation theory. 

Simulations of QCD + QED 

Methods to treat QCD + QED were pioneered some time ago [49]. The total lattice action is 
now a sum of three terms: S = Sym + Sfemtion + Sj7(iV Gauge invariance requires the U(l) field 
to be coupled to fermions in S/ erm(on through Abelian links, [/^ = exp(iqAn). Gauge invariance 
of the photon action Sj/m permits either compact or non-compact photon fields. A non-compact 
formulation of the U(i) gauge field removes photon-photon interactions, so that the theory de- 
scribed by S(7(i) alone is a free theory. As a consequence, photon gauge configurations can be 
easily generated. Non-compact photon fields, however, require gauge fixing. 

Notorious challenges of including QED interactions on a lattice are finite size effects. As 
electric charges are not confined, the resulting long-range interactions suffer from being on a torus. 
For classical electromagnetic fields, field quantization emerges due to the intrinsic properties of the 
torus, namely the flux into any area A in the (z'-j)-plane must be the flux out of the area L,Ly — A, 

7 Naively electromagnetic radiation involving the up-quark alone is 4x more important than that of the down- 
quark. Accounting for photon radiation at the quark level, however, is more involved, see [46], where some surprising 
implications for isospin splittings and Dashen's theorem are explored. 
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because the torus forms a closed surface. Quantized field strengths still lead to finite-size artifacts 
due to the holonomy of the gauge field. These issues with classical electromagnetic fields on a 
torus are mirrored by related issues with their quantum counterparts. Because a torus can leak no 
flux, Gauss's law on the torus becomes trivial. The charge enclosed in the entire lattice volume 
vanishes by periodicity, Q = f^J^dx = Jq V • £ dx = 0. In simulations, this rather extreme finite 
volume effect is mitigated by deleting certain zero modes from the photon field. This non-local 
procedure amounts to forcing the longest-range components of the photon to be non-dynamical. 
We summarize approaches taken so far by groups studying the combined QCD + QED dynamics. 

1) . Since the pioneering work of [49], the first serious attack on the computation of QED 
corrections to QCD has been on domain-wall fermion lattices generated by the RBC collabora- 
tion [50]. In these studies, non-compact QED fields were generated, and domain- wall propagators 
inverted on post-multiplied gauge field configurations, i.e. quenched QED. The QED fields were 
fixed to Feynman gauge with spatial zero modes removed, A^(p = 0,^4) = 0. The physical electric 
charge \e„hy S \ was employed, and a reduction of noise was achieved by averaging over positive and 
negative values of the charge, e = ±e p / !VS . In the most recent study, ensembles with 2 + 1-flavors 
have been used, with pion masses ranging from 250-700 Me V. Results were obtained at a single 
lattice spacing, a = 0.11 fm, and on two volumes V = (1.8 fm) 3 and (2.7 fm) 3 . Finite volume, 
partially quenched chiral perturbation theory + QED was used to account for systematics. 

2) . At previous lattice conferences, reports of QCD + QED calculations using staggered 
lattices have been reported by MILC [51]. These computations use 2 + 1-flavor asqtad con- 
figurations with non-compact photon fields in the quenched QED approximation. The photon 
fields have been fixed to Coulomb gauge, with the following zero modes deleted: A{p^ = 0) = 
and A/[{p = OjXj) = 0. To improve precision, multiple values of the electric charge were used, 
e = ±|e p / 13 , s |, ±2\e p h ys \. Pion masses range from 270-520 Me V, on lattices with volumes (2.4 f m) 3 , 
(2.52 fm) 3 and (3.0 fm) 3 , and lattice spacings, 0.12 fm, 0.09 fm and 0.125 f m, respectively. Con- 
trol over systematics was sought using rooted staggered chiral perturbation theory with QED. 

3) . At the previous lattice conference, reports of QCD + QED simulations were given by the 
BMW collaboration [47], with further updates provided at this conference [53]. These computa- 
tions use an extensive ensemble of 2 + 1-flavor BMW lattices with quenched non-compact QED. 
The photon field is gauged fixed to Coulomb gauge with zero modes deleted; the physical value of 
the electric charge is used. Pion masses range from 135 to 422 Me V, for which an interpolation to 
the physical point can be made, and volumes range from (2.0 f m) 3 to (5.5 f m) 3 . For these compu- 
tations, a single lattice spacing has been utilized so far, a = 0.12 fm. Because of chiral symmetry 
breaking present for Wilson quarks, isospin asymmetric quark masses are radiatively generated at 
€>'{af. s a~ x ). This additive mass renormalization is removed by retuning to the isospin symmet- 
ric point. Finite volume effects appeal - to be well described by a power-law dependence on the 
volume, oc L~ 2 , rather than predicted from chiral perturbation theory + QED. Partially quenched 
chiral perturbation theory is employed to estimate the effect of unquenching QED. While correc- 
tions to individual hadron masses are estimated at the sub-percent level, corrections to meson mass 
splittings are estimated at ~ 5%. 

At this conference: Results of QCD + QED simulations were discussed by three groups. We 
heard of a preliminary study of the electromagnetic corrections to pseudoscalar decay constants 
of light-light mesons and heavy-light mesons [52]. Unique to this study is the use of compact 
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photon fields. In another talk, we heard updates from the BMW collaboration [53], which we have 
included in our summary above. 

We also heard preliminary reports from a comprehensive study by PACS-CS [54]. This study is 
the first to report success at re-weighting configurations to include dynamical QED. Re-weighting 
configurations to include QED effects appears to be suggested first in [55]. Statistical noise in the 
re-weighting factor is shown to be a considerable hinderance — even for a test case of re-weighting 
from e = to e = ^|r. On a QCD ensemble of size 32 3 x 64, the noise was overcome by generating 
QED configurations on a lattice of twice finer resolution, 64 3 x 128, and subsequently averaging 
over all paths in a unit cell to form QED gauge links on the 32 3 x 64 lattice. The masses of the 7t + , 
K + , K°, and D, are used to re-weight PACS-CS configurations to 1 + 1 + 1-flavor QCD + QED 
configurations. The study neglects running of the QED coupling, and QED renormalization of the 
quark masses; however, these are tiny effects. 

4. Outlook 

There have been numerous developments in the study of lattice QCD with classical and quan- 
tum electrodynamics. By the next lattice conference, we might see: further studies of hadron 
correlation functions in weak, uniform external fields; detailed exploration of the QCD phase di- 
agram in strong magnetic fields; computations from multiple groups of the isospin splittings of 
mesons and baryons; and precision determination of the light quark masses accounting for strong 
and electromagnetic isospin breaking. 
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